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The evolution of the curvature perturbation after multi-field inflation is studied in the light of 
the curvaton mechanism. Past numerical studies show that many-field inflation causes significant 
evolution of the curvature perturbation after inflation, which generates significant non-Gaussianity 
at the same time. We reveal the underlying mechanism of the evolution and show that the evolution 
is possible in a typical two-field inflation model. 

PACS numbers: 98.80Cq 



I. INTRODUCTION 

The primordial curvature perturbation Q{k) is strongly 
constrained by observation and provides a unique window 
on the very early universe It is known to have the 
spectrum V^{k) ~(5x 1CT 5 ) 2 with spectral tilt n — l = 
dhxV^/dhik ~ —0.04, and in future one could detect the 
running dn/dlnk as well as non-Gaussianity signaled by 
the bispectrum and trispectrum. 

The process of generating £ begins presumably during 
inflation, when the vacuum fluctuations of one or more 
bosonic fields are converted to classical perturbations. 
Within this general framework, there exist many propos- 
als . 

One proposal is to use two or more inflaton fields, 
which drive inflation in the multi-field model. That 
paradigm has been widely investigated, but it has usu- 
ally been supposed that £(x, t) evaluated at an epoch 
tend just before (or sometimes just after) the end of in- 
flation is to be identified with the observed quantities in 
the spectrum. For this reason, a great deal of effort has 
gone into the calculation of the spectrum, bispectrum 
and trispectrum of C at the end of inflation 043 ■ 

The evolution after many- field inflation has been stud- 
ied numerically in Ref. II Oil using the statistical distri- 
bution of the parameters [ill [l2| . Later in Ref. [l3| the 
evolution of the non-Gaussianity has been investigated. 
In these studies it has been found that there is a minimal 
number of the inflaton field Nf > 10 3 , which is needed to 
realize the late-time creation and the domination of the 
curvature perturbation. Also, the number Nf has been 
related to the creation of the non-Gaussianity. On the 
other hand, the calculation is not analytic and it is not 
clear if the evolution is possible in a two (or a few)-field 
model. 

In this paper, we point out that the actual calculation 
of the curvature perturbation might well depend on the 
evolution after multi-field inflation, even if the number 
Nf is not large. We show that the minimum number is 
Nf = 2, simply because the mechanism requires isocur- 
vaturc perturbation. 



Just for simplicity, consider Nf = 2 with the light 
scalar fields (<f>, a) during inflation. The adiabatic and 
the entropy directions of multi-field inflation are defined 
using those fields. Basically, the "inflaton" (the adiabatic 
field) is not identical to <fi, even if a plays the role of 
the curvaton. The mixing is negligible when a is much 
lighter than <fi; that is the limit where the usual curvaton 
scenario applies. 

Alternatively, it is possible to consider the opposite 
limit, where the fields have nearly equal mass 1 . Can the 
curvaton mechanism work in that limit? A naive spec- 
ulation is that the biased initial condition (o~/4> <C 1) 
might lead to the curvaton mechanism in that limit. In- 
deed the speculation is correct; however to reach the cor- 
rect conclusion we need quantitative calculation of the 
curvaton mechanism in the equal-mass limit. The cal- 
culation details are shown in the Appendix. The usual 
curvaton mechanism is reviewed in SecJIIl and the non- 
linear formalism of the curvaton mechanism is reviewed 
in Sec lIIIl The basic idea of the equal-mass curvaton 
model is shown in Sec lIVI for two-field inflation. Devi- 
ation from the equal-mass limit and the applications are 
discussed in SeclVl 



II. CURVATON MECHANISM 

In this section we review SN formalism used to calcu- 
late £. To define £ one smooths the energy density p on 
a super-horizon scale shorter than any scale of interest. 
Then it satisfies the local energy continuity equation, 



dp{x,t) 
dt 



da(x, t) 



a(x,t) dt 



(j>{x,t)+p(x,t)), (1) 



where t is time along a comoving thread of spacetime and 
a is the local scale factor. Choosing the slicing of uniform 



1 In Ref. 10-13], statistical distribution of the inflaton mass has 
been considered for N-flation. The deviation mMaxMrnin ^$ 
O(10) will be considered in this paper. 
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9, the curvature perturbation is £ = 5 (In a) and 



dt 



= s 



p(t)+p(x,t) 



(2) 



If p is a function purely of p, one will find £ = 0. That 
is the case of single field inflation when no other field 
perturbation is relevant. The inflaton field <fr(x, t) deter- 
mines the future evolution of both p and p. Similarly, the 
component perturbations Q are conserved if they scale 
like matter (p m oc a -3 ) or radiation (p r oc a~ 4 ). 

During nearly exponential inflation, the vacuum fluc- 
tuation of each light scalar field <pi is converted at hori- 
zon exit to a nearly Gaussian classical perturbation with 
spectrum (H/2tt) 2 , where H = a(t)/a(t) in the unper- 
turbed universe. Writing 



C = S[]xi(a(x,t)/a(ti)] = SN, 



(3) 



and taking t* to be an epoch during inflation af- 
ter relevant scales leave the horizon, we define 
N((pi(x, t»), <j>2{%i t*), ■••,£,<*) so that 

COM) = Ni5<f>i(x,U) + -NijSfcfatJSfafaU) H , 

(4) 

where a subscript i denotes d/d4>i evaluated on the un- 
perturbed trajectory. We find 



2y.N i N i r]i i 2 

77 — 1 = ^ I 3 " J — If 

V iV 2 M 2 T N 2 



(5) 



riV = M 2 Vij/V, e = M 2 J2V^/V 2 , (6) 



where M p is the reduced Planck mass. 

The standard curvaton model [3, El assumes that 
these expressions are dominated by the single 'curvaton' 
field ex, which starts to oscillate during radiation dom- 
ination at a time when the component perturbation £ CT 
has negligible contribution to the curvature perturbation. 
Then the non-Gaussianity parameter is given by [l6l . \vf\ 



Jnl 



5 

4r„ 



f 



9 9 

o 

9 



6 



(7) 



where g(a) is the initial amplitude of the oscillation as a 
function of the curvaton field at horizon exit [161]. Here 
r a is identical to ri, which will be defined in this paper 2 . 



III. NON-LINEAR FORMALISM AND THE 
EVOLUTION OF THE PERTURBATION 



curvature perturbation is defined in Ref. [17J, uM as 



dpi 



SN+ I H- 

lp. 3(1 + w i )p l 



SN - 
SN- 



1 



3(l+«;i) \P 
3(1 + Wi) pi 



(8) 



where Wi = 1/3 for the radiation fluid and Wi = for 
the matter fluid. Here a bar is for a homogeneous quan- 
tity, and the curvature perturbation of the total fluid 
should be discriminated from the component curvature 
perturbation Q. The quantity 5pf° = pi — pi in Eq.© 
is the isocurvature perturbation (the fraction perturba- 
tion that satisfies J^Spf = 0), which is defined on the 
uniform density hypersurfaces. 

In order to formulate the evolution of the curva- 
ture perturbation, which is caused by the adiabatic- 
isocurvature mixings, we need to define first the "starting 
point" perturbations at an epoch. 



A. The primordial perturbations 

For the first step, we define the primordial quantities. 
In this paper the quantities at the end of inflation are 
denoted by the subscript "end" , while the corresponding 
scale exited horizon at t*. The subscript "*" is used for 
the quantities at the horizon exit. For our purpose, we 
define the primordial curvature and isocurvature pertur- 
bations at the end of the primordial inflation. 

We find from Eq.©; 

n- - n. e Kl+«>i){Ci-SN) 
~ Pi+3(l + Wi)(f°Pi 

= p> + $p7°- (9) 

Then wc find from p tot = p\ + p 2 = f>\ + pi'- 

J ie 3(l+«'i)(Ci-5A r ) + (1 - f t \ e 3(l+™ 2 )(C 2 -<5A0 = 1; n ) 

where the fraction of the energy density is defined by 

h 



h 



Pi + P2 



(11) 



Expanding Eq. (|10l) and solving the equation for SN, we 
find at first order 1171 



In this paper we consider a clear separation of the adi- 
abatic and the entropy perturbations in a two-field infla- 
tion model. The non-linear formalism for the component 



2 In this paper we use (<f>i,<f>2) for two-field inflation, instead of 
using the conventional (rr, cj>) in the curvaton scenario. 



SN = nCi + (1 - ri)£ 2 

= her + a - nxr] +c 



>adi 



3(1 + wi)pi 



3(1 + wi)pi + 3(1 + w 2 )p2 



(12) 



where Q s ° denotes the second component in Eq.©. r\ is 
defined by 



(13) 
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Defining the primordial adiabatic curvature perturba- 
tion (C mf ) just at the end of inflation, the component 
curvature perturbation (Q) can be split into £ mf and Q so . 

The obvious identity is 



?"l,endC^ond + (1 - r l,cnd)C2^nd — 0; 



(14) 



which is valid at the end of inflation. Apart from that 
point the deviation due to the evolution of r\ becomes 
significant. 

The parameter of the fluid (wi) is constant when pi 
behaves like matter (wi = 0) or radiation (wi — 1/3), 
and a jump (e.g, Wi — — > Wi — 1/3) is possible when 
instant transition is assumed. In this paper we are using 
the sudden-decay approximation for the curvaton mecha- 
nism. 3 We also assume that the inflatons start sinusoidal 
oscillations just at the end of slow-roll. 

The curvature perturbation in the standard curvaton 
scenario is usually expressed as 



5N = nCi + (1- r!)C inf . 



(15) 



Assuming that Cf° > C > Cf °j one will find Ci ^ Cf° 
and C2 ^ C inf , which gives Eq.^ from Eq.(H2|). Usually 
the above approximation is justified when mi < m-2 and 
the curvaton is negligible during inflation. 

In this paper we are considering the equal-mass limit 
(mi ~ m 2 ), which is in the opposite limit of the conven- 
tional curvaton. In Appendix we show the validity of the 
above approximations and derive the quantitative bound 
on the ratio 4>i/4>2- 



IV. A BASIC MODEL 

In this section we show why the curvaton mechanism 
can create the dominant part of the curvature perturba- 
tion after conventional chaotic multi-field inflation, nei- 
ther by adding extra light field (curvaton) nor by intro- 
ducing many inflatons. The calculation clearly explains 
why and how the curvaton mechanism works in the equal- 
mass limit (mi ~ mi). 

We assume (for simplicity) that after inflation the field 
4>2 decays immediately into radiation and 4>\ starts sinu- 
soidal oscillation at the same time. Then cf>i decays late 
at Hdi <C Hi. There is no mixing between these com- 
ponents. Here Hi denotes the Hubble parameter during 
primordial inflation. 

In this scenario, we consider two phases (^4, B) charac- 
terized by wia = and wib — 1/3. Here the subscripts 
A and B denotes the quantities in the phase A and B. 
They are separated by the uniform density hypersurface 
Hdi — T\\ 



3 Authors of Ref. Il9|| consistently accounted the curvaton decay 
when the curvaton decays pertubatively and showed that the 
correction to the potential can be significant. 



(A) pi; oscillation, pi] radiation 
(wi = 0, w 2 = 1/3) 

(B) Radiation 

(wi = w 2 — 1/3). 

The important assumption of the model is that the 
transition occurs on the uniform density hypersurfaces 
so that we can neglect additional creation of 5N (modu- 
lation) at the transition. 

We find in phase (A); 



6N = tiaCia + (1 - ruOCa, 



(16) 



where the subscript "A" (or "i?") is omitted for (2, since 
£2 is constant during the evolution. Here we used the 
definition 



riA 



3pi 



3pi + 4p 2 



(17) 



Consider a simple double-quadratic chaotic inflation 
model in the equal-mass limit. The potential is given by 



\m- (4 



) = \™ 2 &, (18) 



where ^1,2 are real scalar fields. Besides the potential, we 
need the interaction that causes difference in the decay 
rates. Fig[T] shows the evolution of the densities after 
inflation. The end of chaotic inflation is given by 



^l,end 



^2, end — ^r,end 



Mi 



(19) 

Since the potential is quadratic during inflation, we find 

1 <5</v* 



inf 



(20) 



In this section we consider 9 <C 1 , which leads to the sim- 
plifications sin# ~ 9 and cos 9 ~ 1. Our approximations 
are based on the exact calculation in Appendix A. 




FIG. 1: tend, td2, iosc and tdi denote the time at the end 
of inflation, 02 decay, the beginning of 4>i oscillation and <f>i 
decay, respectively. Our scenario is shown in the left-hand 
side, which gives the time-ordering t en d — t osc < tjs < tdi- 
The usual curvaton scenario is shown in the right-hand side, 
which gives t end < t d 2 < tosc < tdi- 
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1 2 



Using Eq. ([25)l . we find 



— f- < 5i] x 10~ 3 , 

Mr, 



(27) 



which does not always require significant suppression. 
The ratio y = y/T 1 /r 2 is calculated in Eq. ([A24l) and 
is given by 



V ■ 



\fNLt 



(28) 



FIG. 2: The straight dotted line with an arrow is the inflaton 
trajectory, and the circle gives the uniform-density surface 
along which the entropy perturbation 8s appears. 



From Eq. (|A19[) , we find the component perturbation of 
the late-decaying component {<j)\ ) at the end of inflation: 



c 



1A 



i w 



end 



3 Pi, 



end 



2 59 

3~T 



(21) 



The usual approximation of the curvaton mechanism is 
Ci% 3> C ■ The validity of this approximation is exam- 
ined in the Appendix. 

From Eq. (|A2ip , the final curvature perturbation is 



59 1 fS9\ 

T + 2 T 



(22) 



Defining the ratio y = y/Ti/T2, r±- (r\ evaluated in the 
phase (A) just before the decay) is given by 



39 2 



39 2 + Ay 



(23) 



The non-Gaussianity parameter has been calculated in 
Ref. [H. We find for 9 < 1: 



5 

4ri 
5 



g"g 



f NL ~ — li + z-f 



r-i 

3 6 



4ri_ 



(24) 



Further simplification is possible when 59 = <5s*/(/> r * 
and Vss, = T-'8<t> r , ■ For the quadratic potential we have 
<p r * = 2y/N e M p 3> e ,ond, where N e is the number of 
c-foldings during the primordial inflation spent after the 
corresponding scale exits horizon. The condition of the 
curvaton mechanism £ fin > ( lnl gives 



< 3 n - V ' 



5 1] 



6 fNL 



(25) 



Here 9 should be less than 1 but does not require many 
orders of magnitude. From the CMB spectrum we find 
the normalization given by 



V 



1/2 



n- 



Hr 



6ttJW:9 m v 



5 x 10" 



(26) 



We thus find that the difference between <f>\ and <\>i decay 
rates is in the conceivable range. 

The above conditions tell us how small 9 and y have to 
be to get a given CMB spectrum and Jnl- They have to 
be some orders of magnitude below 1 but not very many. 

If the potential during inflation is both symmetric and 
quadratic, we find 77 = 771 =772- We thus find the spectral 
index 



1 = -2e + r/ = 0, 



(29) 



which shows that the above model requires deviation 
from the symmetric potential. 

Lookingback into the many-field inflation, the model 
in Ref. [T(| assumed that the inflaton masses are not 
exactly the same but may have statistical distribution 
around the mean value. In that case, the cancellation 
in the spectral index is not realistic. Since the deviation 
from the symmetric potential is expected, we need to 
examine what deviation is needed for the model. Then 
we can understand why and how the curvaton mechanism 
works in the many-field inflation model. 



V. DEVIATION FROM THE SYMMETRIC 
POTENTIAL 



The deviation from the symmetric quadratic potential 
can be classified as follows; 

1. A small mass difference (1 < 1712/1711 ^10) 

The spectral index does not vanish when the double 
quadratic potential has different (but not so much 
different as the usual curvaton) mass. The slow-roll 
parameters are 



H 2 



3Hf 



(30) 



where the fraction of the density is given by fi = 
. The spectral index is shifted from n s — 1 = 
and is given by 

n s — 1 = -2e H + 2i]i 

~ -2[7 ?1 / 1 + 7 ?2 (l-/i)]+2r ? i 



5 



-2(ry 2 -77i) 
-2P V2 
P 



(31) 



where P : 



< 1. The observation 120] shows 



n s -l = 0.037± 0.014, which suggests N e < 40 and 
requires secondary inflation (2lj . 

Besides the spectral index, mi < m2 suggests that 
the oscillation of the field 4>i is slightly delayed com- 
pared to (f>2- The delay may enhance the density of 
(f>i at the beginning of the oscillation, while the ini- 
tial pi density may be reduced since mj is smaller. 
Defining y c g = y/Ti/m± and 6 c g = <f>x/<t>2 at the 
end of inflation, we find 



riA- 



3mf 



2Aj 



3mf^ + 4ml0l(g) 



2/cff 



36l off 



(32) 



which gives a similar bound for 9 c s (y c s) 



2. Heavy curvaton (mi > 7712) 

Usually the curvaton is assumed to be much lighter 
than the inflaton; however this assumption could 
be avoided. We consider the curvaton mechanism 
when the curvaton is slightly heavier than the in- 
flaton. 

We assume pi > p\. Once it is assumed at the 
beginning of inflation, it remains true during infla- 
tion. 4 Then </>i-oscillation starts during inflation. 
It begins when 



m\ = Hi 



6M 2 



(33) 



where the subscript "osc" denotes the beginning 
of 0i-oscillation. From the above equation and 
^2 1 osc — 2y/N2M p , where N2 is the remaining 
number of e-foldings after the beginning of 0i- 
oscillation, we find 



No, 



3m 2 
2mf ' 



(34) 



= e 3 ^ 2 
we can estimate 



Defining y cS = e iN * ^/^l7^2 and 



[0i M 



riA- 



2 

osc 



3^2 SC + 42/eff 



2 J osc ) 



(35) 



Unfortunately, the spectral index is 

n s — 1 ~ ~2en + 2r\\ 
~ -2?72 + 2r/i 
~ 2771 > 0. 

3. Symmetric but Non-quadratic 



(36) 



The potential could be dominated by a polynomial 
V(cj> r ) oc </>p at the moment when the perturbation 
exits horizon, while it can be approximated by the 
quadratic potential during the oscillation. For the 
polynomial we find <f> r * ~ \J2pN e M p and the slow- 
roll parameters 



2 p e 



M, 



;P(P- 1) 



(37) 
(38) 



The spectral index is shifted and is given by 



M, 



<t>l 

p-2 

2N r< 



f \p 2 ~2p(p~l) 



(39) 



The result suggests that p < 2 is needed for the 
scenario. In that case the mass and the coeffi- 
cient of t he p olynomial must run in the trans- 
Planckian [22]. p — 1 would correspond to mon- 
odromy in the string theory and it requires N e < 
20. p < 1 is an interesting possibility if the effective 
action allows fractional power. 



A. A model with a complex scalar 

An inderesting application of the idea is that a conven- 
tional 2-field multiplet contains both inflation and the 
curvaton at the same time. Consider a complex scalar 
field $ = 02 + *0i, which gives the symmetric potential 

y ( $) = I m 2 |<i> 2 | 2 = I TO 2 ( 2 + 2 ). (40) 

First, consider a small symmetry breaking caused by 
A 4 / $ + $* x 2 



AV 



M 2 



(41) 



where A <C M is assumed. </>2-oscillation may cause sig- 
nificant particle production when there is the interaction 
given by 



Ant = g{® + $*)#, 



(42) 



4 The opposite condition (pi > P2) requires 0i» > M p , which 
suppresses the component perturbation of the curvaton and does 
not realize the curvaton mechanism. 



which can lead to significant ^-production at the en- 
hanced symmetric point (02 ~ 0) [23| . The coefficient 
of the interaction could be small (g ~ A/M <C 1) when 



6 



it is suppressed by a cut-off scaie. ip may decay quickly 
into radiation since the amplitude of the oscillation after 
chaotic inflation is very large [23| . 

Define 8m? = ^p-. If 8m? is much smaller than m 2 , 
the cancellation in Eq. ([3!l is still significant. On the 
other hand, it is possible to assume 8m? I m? ~ 0(1) 
(which is still within the conventional set-up of multi-field 
inflation) one obtains P ~ 1 and n s — 1 ~ — -A-. Again, 
the scenario requires additional inflation stage |2l| . 

Second, consider the case in which the potential dur- 
ing inflation is dominated by a polynomial V($) oc 4> p . 
The curvaton can dominate the spectrum, however the 
spectral index becomes 



_Ml 

p-2 
2Y e ' 



[p 2 - 2p(p - 1) 



(43) 



The scenario requires p < 2. 



B. Sneutrino inflation 

It is possible to assume small inflation "before" the 
multi-field inflation. The observed spectrum of the curva- 
ton perturbation exits horizon during the first inflation. 
In that case ch is determined by the first inflation and 
the cancellation in the spectral index is avoided. This 
scenario uses multi-field inflation for the curvaton infla- 
tion m. 

The usual sneutrino inflation [27] uses m ~ 10 13 GeV 
to satisfy the CMB normalization. When the condition 
is combined with the gravitino problem, Yukawa cou- 
pling of the first generation sneutrino (single-field infla- 
ton) must satisfy (iy^v)n < 10~ 12 , whilst other Yukawa 
couplings will not be so small. Here Y v is the neutrino 
Yukawa matrix. 

In this section we consider multi-stage inflation, in 
which three sneutrinos play crucial role. We assume that 
the first single-field inflation is caused by the third gen- 
eration sneutrino, and the secondary two-field inflation is 
caused by the first and the second generation sneutrinos 
with the mass M\ = M 2 = M. We assume M 3 > M for 
the third generation. 

The reheating after two-field inflation is due to the 
decay of the second generation sneutrino, which gives 
the reheating temperature 



T R = 



90 



1/4 



where the decay rate is 



Ti ~ — (Y v Y!j)..M. 



(44) 



(45) 



From Eq. (|A9[) . the curvaton mechanism is significant 
when T 2 3> Y\. For the two-field sneutrino inflation, 



which is the secondary inflation of the above scenario, 
we find 

Cl < **{ ( y v y2 )ii ) m p (46) 

Here the mass of the first (second) neutrino is 

< H u > 2 



{m v )ii ~ (Y v Yj)i 



M 



(47) 



We thus find for the given neutrino mass (m„)n and 
{m v )22\ 



v i/2 If {m u ) 22 \ 
^ 37rV(m„)iJ 



1/4 



M 



(48) 



The reheating temperature after inflation is given by 



T R = 



( 45 



1/4 



M 



: ^/(m v ) 22 Mp, (49) 



while the temperature just after the curvaton decay is 



_ f 45 \ 1/4 VI 



-^(m u ) n M p . (50) 



R \8w*gJ <H U > 
We may write the spectrum using and T' R ; 

1/2< J_ (TnV /2 ^ 
Cl 3tt \T' r J M p 



(51) 



When the primary inflation gives the number of e- 
foldings Ni, the spectral index is 



n, — 1 ~ — 2ejj — • 

Ni 



(52) 



The observation gives n s — 1 = 0.037 ± 0.014, which sug- 
gests 20 < N x < 40 for the first inflation. 



C. N-flation 

The two-field inflation model considered in this paper 
is a simplification of the N-flation model [24j. The In- 
flation has been studied using statistical argument [To| . 
which helps us understand the results obtained above for 
the two-field model. 

Assuming (for simplicity) the same potential for all Nf 
fields, we find 



Using the adiabatic field defined by 

N f 



(53) 



(54) 
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we find the potential 



VI. CONCLUSIONS 



V(4> r ) = -m 2 4> 2 r . 



(55) 



If we assume uniform initial condition <j> n ~ <pa, the model 
is identical to the two-field model with 9 ~ 1/ \fWf <C 1. 

For the number of e-foldings N e ~ 60, the usual cur- 
vature perturbation created at the horizon exit is given 
by 



C inf = _ Hl 5 4L 



= 2N B 



(56) 



where Hf = Nf 6 ™ff° is the Hubble parameter during the 

primordial N-flation. 

Suppose that the decay rate T n is uniform except for 
a field (j>i, which has Ti <C T n . Here the density ratio 
becomes r* ~ Repeating the same calculation, we 
find 

Ci ^^Pi = 2^i ^ 2 

3pi 3 00 3 4> r 

V 1 ^ < V^ 2 is possible when N f > iV e 2 . This gives 
the minimum number of the fields that is needed for the 
curvaton mechanism and it explains the numerical calcu- 
lation in Ref. (loj . 

In the above scenario, the curvaton is one of the in- 
flaton fields that are equally participating l/Nf of the 
inflaton dynamics. 

At the end of inflation, the fraction of p\ is 

ri{t end ) = -jr- < 1, (58) 
N f 

while at the decay of <f>\ it can grow; 

ri(t d ecay) = I (tend) x 



1/2 



(59) 



We need for the curvaton mechanism (i.e, (^-domination) 



3 V S (j> r Nf 



1/2 



> 2N, 



TV 



(60) 



which leads to 



r„ 



1/2 



< 



3A^ e 



(61) 



Significant non-Gaussianity (/jvx) requires r(td eC ay) 
0.1, which gives 



1/2 



10 

w 



(62) 



If the distribution is statistical for the decay rate, we 
need Nf ^> 1 for the strong suppression (Fi/T n -C 1). 

In this section we found that the evolution after infla- 
tion may dominate the curvature perturbation when Nf 
is large. Our result explains the numerical calculation in 
Ref. 53 . 



The evolution after multi-field inflation can change the 
curvature perturbation. In this paper we considered a 
conventional two-field inflation model and showed that 
the curvaton mechanism after multi-field inflation could 
be significant when the decay rates are not identical 5 . 
Interestingly, the mechanism works for a complex scalar 
field $ = <f)2 + i<f>i. 

The previous numerical study [13 showed that Nf 3> 1 
causes significant evolution of the curvature perturbation 
after inflation as well as the creation of significant non- 
Gaussianity. We showed that the same is true for two- 
field inflation, in which 9 <C 1 is required instead of Nf 3> 
1. 

The source of the curvaton mechanism is the entropy 
perturbation generated during multi-field inflation. Since 
the uniform density surface of the multi-field potential 
is fiat by definition, the perturbation on that surface is 
inevitable. 

Our results suggest that many-field inflation must be 
considered with care. A large number (Nf > 10 3 ) can 
easily explain the required condition for the curvaton 
domination. 
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Appendix A: Calculation details 
1. Evolution of the curvature perturbation 

In this Appendix we show the calculation details of the 
evolution after inflation. 

We first assume that the potential is quadratic and 
symmetric during chaotic inflation. In our formalism £ 
is defined at the end of inflation. The entropy perturba- 
tion is realized by 89, which is the perturbation of the 
angle 9 in FigJ5J 

The spectrum of the entropy perturbation during in- 
flation is Vgs, — (Ht,/2n) 2 . The entropy perturbation 
causes the fraction perturbation between densities. Using 
69, the densities of the components and the isocurvature 



A similar but another story has been discussed in Ref. 128 
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perturbations at the end of inflation are given by 
1 



Pi, end 



2uend|2 sin 2 ( 



-m 2 |0! 



2 " r 1 ~ 2 

,2 i\/r2 



m 2 M 2 

sin 2 (Af ) 



8pf° nd ~ m 2 Af 2 (sin cos 0)<50, (A2) 

' " " " ' (A3) 

(A4) 



/02,cnd = -m 2 \(j) C r nd \ 2 COS 2 V ~ - 
<K°end + ^fend = 0- 

We find at the end of inflation: 

P 1 ---.2, 



m 2 M 2 „ 
— =— ^ cos 2 I 



,fl 
5/1 



Pi + P2 

50 
= [sin 20 



sin 2 0, 



89 = 2 [sin cos 01 1 



(A5) 



(A6) 



The expansion with respect to 89 makes no sense when 
89/ sin0 > 1 or 89/ cos0 > 1 [29[. We are excluding those 
regions. 

Creation of the curvature perturbation after inflation 
requires the decay rate T\ <C ^2- In the phase (A) we 
find 



C 



2 cos 



I A 



^2A 



3 sin0 
1 sin0 



2cos0 



.60. 



(A7) 
(A8) 



Using Ea.(fl"2"|). we find 

-fin 



2 cos 1 sin 



3 sin0 2 cos 
4ri_cos 2 0-3(l-r-i_) sin ° 



inf 



6 sin cos 



inf 



(A9) 



where rj_ denotes the value of t\a evaluated just before 
the end of the phase (A) . 
The evolution is 



Domination by the curvaton density (ri_ ~ I) requires 
\/T\JT^ < tan 2 0. 

The CMB spectrum requires V ( nn ~ (5 x lCT 5 ) 2 [2(ij . 
The requirement is trivial when £ fin ~ £ mf , 6 while in the 
opposite case C fin > C mf ; m which the curvaton mecha- 
nism dominates, we need the condition 



2 cos 9 1 sin 
— n_ = (I — ri_J = 

3 sin0 2 cos0 



> 



?7</>r* 



(AI3) 



Solving Eq. (|A13j ) for ri_ and using Eq. (|A12| ). we find 

(A14) 



/Ti 2??tan0 - 3tan 2 

— < — = < f. 

T 2 4 + 277 tan 



This equation also shows that 2r) — 3tan0 > 0, which 
gives 



tan0 < -rj. 
3 

The CMB observation gives the normalization 



(A15) 



2 cos 9 1 sin 

3 ri ~ sm~0 ~ 2*- 1 ~ 7 * 1 "' ) cosl 



,1/2 



T^ 2 - 5 x f 0" 5 (A16) 



Defining k = 5x fq-s and y = yjTx/Yz, we can solve 
Eq. dAf 6P for y and find 



2fc - 3tan0 
V ~ 2fc + 4tan _1 
fc0 

~ Y' 

To avoid y < 0, we need the condition 

3 n , 

— tan0 < k. 



(A17) 



(A18) 



The perturbations can be expanded up to second order. 
We find 



Pi- 



92- = 



m 2 M 2 



p sin 2 I 



m 2 M? , 
cos 2 I 



which leads to the ratio 

P2- 

Pi- 



cos 



Odl 

Qd2 
("end 



a<i2 
adi 



— ) (MO) 



(AH) 



Therefore, in the radiation dominated Universe we find 
3pi_ 



n_ = 



3pi_ + 4p 2 - 

3 sin 2 



3sin 2 + 4cos 2 v / rT7T2' 



(A12) 



COS0 



sm( 



ISO 

2 



sin 



89 
cos 



- cos 20 
2 \ sin ( 



- cos 20 , 

2 V cos 



(AI9) 
iA20) 



Using Eq. ([T6|) . the final curvature perturbation after the 
decay is 



fin 



2 cos 1. sin ( 

n-^^ - o(i - ri_)- 



sm f 



cos( 



6 Note however the non-Gaussianity is not trivial because the cur- 
vaton perturbation may still dominate the second-order pertur- 
bation 13111 . 
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1 cos 29 1. ,cos26» 



3 sin 

inf 



cos 2 9 



{sef 



< 

4ri_ cos 2 6 - 3(1 - n_) sin 2 



cos 26* 
+ 3sin 2 20 

+c inf 

4ri_ cos 2 ( 



3 sin 29 

[4n_ cos 2 - 3(1 - n_) sin 2 6] (<5<9) 2 
-3(l-ri_)sin 2 

(A21) 



3 sin 26* 



cos26> 



sin26» 



^inf 



When the curvaton perturbation dominates (9 <C 1), 
the non-Gaussianity of the spectrum is measured by 



5 cos 29 



NL 



4ri_ cos 2 6 - 3(1 - r x _) sin 2 



(A22) 



Using Ea. (|A12J) . we can substitute r X - in Eq. (|A22D . 
Then solving the equation for y, we find 



- tan 



4 cos 2 9 

5 cos 20 



(A23) 



Barring cancellation, the above equation gives a simpli- 
fied formula 



5 



(A24) 



Being combined with Eq. (|Al7j) . which has been obtained 
using the CMB normalization, we find 



NL C 



(A25) 



We thus find (from /jv l and CMB using the definition of 

k) 



v l/2 



6 x 1(T 5 x f NL 



or equivalently 

Hi ~ 6 x 1(T 3 x f NL 9M p . 
Solving the equation for 9, it gives 



1 



x 10'' 



(A26) 



(A27) 



(A28) 



Using Hj in Eq. (|A27|) and calculating the tensor to 
scalar ratio r g , we find [301 ] 



f NL 9 2 x 10 4 . 



(A29) 



Considering the natural bound T2 < Hi and Ti > 
H nU c, where -ff nuc is the Hubble parameter at the time of 
the nucleosynthesis, Eq. (|A23[ ) gives the lower bound for 



> 



Hi 



1/4 



(A30) 



Besides the above condition, we have another condition 
coming from 9 > 69. Since we are assuming quadratic 
potential in the trans-Planckian, we have 69 — 5s/<f> r * 

leads to 



and </> r * = 2s/N e M p . Then 9 > 



9 > 0.01 



El 

Mr, 



(A31) 
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